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Abstract 

Numeric model of the bremsstrahlung response of homogeneous mag- 
netoactive plasma on a gravitational wave with e+ polarization was con- 
structed. Electromagnetic response dependencies on the plasma and grav- 
itational wave parameters were determined. 

1 Introduction 

The equations of the relativistic magnetohydrodynamics (RMHD) of a mag- 
netoactive plasma in a gravitational field were formulated in paper using 
the equality requirements for dynamic velocities of plasma and electromagnetic 
fielcS. These equations were obtained on the basis of the Einstein and Maxwell 
equations. Also the remarkable class of exact solutions of obtained RMHD 
equations was found. It explains movement of a magnetoactive locally isotropic 
plasma in a field of plane gravitational wave (PGW). This class was called gra- 
vimagnetic shock waves (GMSW). It describes essentially nonlinear processes 
which do not exist in a linear approximation of the magnetohydrodynamics and 
essentially relativistic processes in terms of domination of massless electromag- 
netic component in a magnetoactive plasma. 

It was shown in paper ^ that the GMSW in pulsars magnetospheres may 
be the highly effective detectors of gravitational waves from neutron stars. Par- 
ticularly, giant pulses which sporadically appear in radiation of some pulsars 
may be the observation result of transferring energy from a gravitational wave 
to the GMSW. Estimations made in [5], [5] allow to connect giant pulses in 
radiation of pulsar B0531-I-21 with gravitational radiation in the basic mode 
of oscillations from this pulsar. In fact, at this moment it is sufficiently diffi- 
cult to say about identification of giant pulses as the electromagnetic display 
of gravimagnetic shock wave evolution in pulsar magnetosphere and to connect 
unambiguously these pulses with pulsars gravitational radiation. Nevertheless, 
the idea of analyzing the influence of gravitational waves from a compact astro- 
physical object on its own electromagnetic radiation is highly productive for 
solving the problems of gravitational waves detection. 

In fact, the main difficulties of gravitational waves detection in the Earth 
conditions are: 



^Till 2000 Yu.G. Ignatyev wrote his name as Yu.G. Ignat'ev. 

■^This requirement is completely equivalent to the condition of plasma infinite conductivity, 
see Ref. [l]. 
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1. Negligible amplitude of gravitational waves in the Earth conditions {h < 
10""'^^) because of significant distance between relativistic astrophysical 
objects and the Earth. 

2. Sporadic nature of events leading to radiation of gravitational waves inside 
relativistic astrophysical objects with enough power. This does not allow 
to connect unambiguously received signal with a fact of gravitation radia- 
tion detection. 

3. Impossibility to construct relativistic detectors with anomalous highly ef- 
fective parameters for registration of gravitational waves in conditions 
of the Earth laboratory (super-strong magnetic fields, highly anisotropic 
working body of detector, low level of background noise, etc.). 

It is possible to avoid these problems if one can transfer detector directly 
close to a relativistic astrophysical object. In this case one always has prepared 
electromagnetic signal and there is no need to convert it in other forms so it 
allows to do correlation analysis. If the detector's working body is the magne- 
tosphere of relativistic astrophysical object, the optimal for gravitational waves 
registration parameters of detector's working body will be achieved automat- 
ically: super-strong magnetic fields, ultrarelativistic equation of state, highly 
anisotropy, etc. 

The fundamental importance of the GMSW for theory as the direct conver- 
sion effect of gravity- waves energy into electromagnetic energy leads to neces- 
sity of more detailed and comprehensive researches. In Ref.[4] the strict proof 
of the GMSW hydrodynamic theory based on relativistic kinetic theory was 
given. In [1], [2], was shown that GMSW realizes in essentially collisionless 
nonequilibrium plasma within anomalous strong magnetic fields. Isotropy of 
a local plasma electron distribution essentially violates in the such conditions 
due to strong bremsstrahlung. Therefore, an anisotropy factor of magnetoac- 
tive plasma is highly essential for effectiveness of GMSW formation mechanism. 
Hydrodynamic model of GMSW in anisotropic plasma with adjusted correla- 
tion between parallel and perpendicular components of plasma pressure was 
constructed in [S]. It was based on the general equations of RMHD. Partic- 
ularly, in [5] was considered the elementary linear correlation. The research 
made in [5] discovered the strong dependence of GMSW effect upon the plasma 
anisotropy degree. That fact led to necessity of constructing the dynamic model 
of anisotropic magnetoactive plasma movement in a gravitational radiation field. 

Process of a gravitational wave energy pumping over into the electromagnetic 
energy is describing with the help of the energobalance equation introduced in 
Ref.[T]-[5]. It performs the fact of total momentum conservation law inside the 
system "gravitational wave -I- magnetoactive plasma" . In paper |6j analytic 
research of essentially nonlinear equation was done and some features of the 
solution were detected. However, because of software existed in 1998 and other 
reasons the total research of GMSW evolution was not done and the parameters 
of bremsstrahlung response of magnetoactive plasma on gravitational wave were 
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not obtained. This paper is dedicated to these problems solution. Here we set 
the unit system where c = G = h = 1. 



2 Gravimagnetic shock waves 

Let us reproduce the main results of the GMSW theory which are necessary 
for the goals of this paper. Let us set the metrics of vacuum PGW with e+ 
polarizatiorQ which propagates alone the Ox^ axis: 

ds^ = 2dudv - L^[e^^{dx^)^ + e-^^{dx')^], (1) 

where /3(u) is an arbitrary function (the PGW amplitude); the function L{u) 
(the PGW background factor) obeys an ordinary second order differential equa- 
tioE0 ; It = -^{t — x^) is the retarded time and v = '^{t + is the advanced 
time. Let in the absence of PGW {u < 0) be given a homogeneous magnetic 
field directed along the Ox^ axi^: 

H,{u < 0) = S^Ho. (2) 

In general, alternating electromagnetic field which appears in the presence of 
gravitational wave has only spacelike magnetic component in the comoving 
frame of reference moving with local velocity [1]: 

= v'' Fki\ Ffez= ^^zimi^'"; (3) 

{Fik - Maxwell tensor, Fki - dual Maxwell tensor, r^kUm - discriminant ten- 
sor). The electric component of electromagnetic field in the comoving frame of 
reference equals zero [Tl: 

= v'^Fk^ = 0, (4) 
therefore the energy-momentum tensor (EMT) of electromagnetic field is: 

f ^ {2H^v,Vj - 2H,Hj - g^jH^) . (5) 

Squared magnetic field strength is determined as [1]: 

H^ = -iH,H) = ^Fi„,F'"'. (6) 

Thus, the trace of the EMT of electromagnetic field is equal to zero: 

T= g'' T^J^ 0. (7) 



•^The case of two polarization states will be considered in next paper. 
*See, for example, Ref. [71. 

^The general case of a magnetic field arbitrarily directed in the plane x^Ox^ was considered 
in the cited papers. 
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Invariants of electromagnetic field comply with the conditions: 



F^kF''' =F''' F^k^ 2H^ > 0; F^k F''' = 0. 



(8) 



Let further magnetoactive plasma be homogeneous but anisotropic in general 
in the absence of PGW. In gravitational field the EMT of anisotropic magneto- 
active plasma is [5]: 



T'^^ {e + p±) v^v^ - p±g'^ + {p\\ - p±) h'h^ 



(9) 



where /i* = W/H - spacelike unitary vector of a magnetic field {{h, h) = — 1), 
P± P\\ - perpendicular and parallel components of the plasma's pressure, and 
according to ([3]): 

iv,h) = 0, (10) 



And the EMT Q according to virial law complies with the condition: 

p 



(11) 



T= e- p|| - 2p_L > ^ j3|| + 2pj_ < e. 
Let's further suppose a barotropic equation of state: 

P\\ = k\\£; P± = k±e, (12) 
where coefficients of baratrops /c|| , k± by reason of (fTT|) are follow the inequality: 

fc||+2/c_L<l- (13) 
Then in a presence of PGW the exact solution of RMHD equations is [5]: 



1 

71 



V2 = 0; Vu^ TT^; 

2v,, 



^2^fe|l+fcig/3(fc||-fej_) 





£ —e 



^l + fcj_2^2(l+fc||)g2/3(fc||-fcj_) 



H = Ho 



AL(i+'=ii)e-'3(i-'^-||) 



-a± 



where 
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G [1,2], 



1 - kj_ 

A{u) is the governing function of GMSW. 



(14) 
(15) 

(16) 
(17) 

(18) 
(19) 

(20) 
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n - local charged particle density, - dimensionless parameter: 



47r(e + 

Variables from above marked with zero are given in the absence of PGW. 
The RDMD equations solution consists of the physical singularity on the hy- 
persurface : u = : 

AK)= [l-a2(e2'5("*)-l)l, (22) 

on which the densities of the plasma energy and of the magnetic field tend to 
infinity, the dynamic velocity of the plasma as a whole tends to the velocity of 
light in the PGW propagation direction. In this case the ratio of the magnetic 
field energy density to the plasma energy tends to infinity. The above singularity 
is the GMSW spreading in the PGW propagation direction at a subluminal 
velocity. According to Eg. ip^ the conditions of the singularity arising are 

/3(u) > 0; (23) 

o? >\. (24) 

The extremely important fact is that, the singular condition is even possible 
in a weak PGW <C 1) on the condition of a highly magnetized plasma 
{o? 1). In this case the singular condition, according to (f^. arises on the 
hypersurfaces u — u.^,: 

1 

It follows from ([Ti| - p7|) that, by /3 > the plasma moves in the GW 
propagation direction (xr" ~ -^(v^ — Vv") > 0), hy /3 < - in the opposite 
direction. 

The singularity was removed by taking into account back influence of the 
magnetoactive plasma on a PGW. It leads to effective absorption of a PGW 
energy by the plasma and to PGW amplitude restriction. The simple model 
of energy balance which describes that process was constructed in p . Gravi- 
tational wave with metrics ([T]) in WKB-assumptior|l corresponds to the EMT 
with one nonzero component: 

Tuu=^{P'?- (26) 

Let /3*(u) be a PGW vacuum amplitude and (3{u) be a PGW amplitude in 
consideration of absorption in plasma. In this case the energobalance equation 
in the short-wave approximation becomes: 



PM = T^- (25) 



°In this case it corresponds to the condition oi p \ = c/uj where p is a space-time 
curvature radius. 
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(/?:)' = + 47r(r„„-T^.n), (27) 

where Tik is total plasma EMT. Under condition of ^ 1 Eg. p^ may be 
written in the form: 

ql^q' + eV{q), (28) 

where q — /3//3o and the dot signifies a derivative in the dimensionless time 
variable s: 

s = V2uju, (29) 

{uj - the PGW frequency), V{q) - potential function which in a weak PGW 
becomes: 

Viq)^A-^a^iq)-l, (30) 
where is so-called the first parameter of GMSW ^: 

- S- 

Eq.(|28|) may be treated as an equation with respect to the variable q. On the 
other hand, (|28p completely coincides in its form with the energy conservation 
law of a 1-dimensional mechanical system described by the canonical variables 
{q{s), q{s)} jOl, where V{q) is the potential, q"^ is its kinetic energy and = 
is its total energy. 

Let us introduce the new dimensionless parameter: 

T = 2a^Pa (32) 
- {the second GMSW parameter) and rewrite (|22l) in a weak PGW as: 

A{q{s)) = 1 - 2a^l3oqis) - 1 - Tq{s). (33) 
It leads from ((55|l : 

<i = -^. (34) 

To analyze the system behavior, let us suppose that the moment s = 
corresponds to the front edge of a GW, whil^ll: 

/3* « /3o(l - cos(s)) ^q^Kil- cos(s). (35) 

According to the system starts with negative value of the governing 

function derivative and with function value equal to 1: 

A(s) ~ — Tsins « — Ts; 

(s^+0). (36) 

A(s) w 1 - T(l - coss) w 1 - T^; 



^This provides zero PGW metrics derivatives at the moment s = 0, i.e. class of metric 
functions. 
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The energobalance equation (1^ according to ([50]) . ([M)) . ([55]) becomes: 



^-4gi - 1 



T2sin2(s). (37) 



Solving the Eg. ([57)) with respect to A we obtain: 



A = T^^sm\s) - e [A-49^ (s) - ij . (38) 

Integrating according to first of all it's necessary to take negative branch 
of the Ea. (j38p but when we reach the minimum value of the governing function 
we should change it by the positive one. From (j38|) we obtain the minimum 
value of the governing function which is reached by s = 7r/2: 

1 

e 



Am^n = ( + 1 ) , (39) 



where: 

The maximum accessible density of a magnetic energy is 



and it does not depend on a plasma equation of state ([T^ . Also plasma ve- 
locity in the GMSW does not depend on equation of state. And the maximum 
plasma energy density without magnetic field depends on the exponent of plasma 
anisotropy: 

/ 1 ^ 3(1+'=^) 



It is maximum for the ultrarelativistic plasma with zero valuation of the 
parallel pressure. 

Thus, the maximum value of the local response amplitude of a highly mag- 
netized plasma (a^ 3> 1) with linear state equations does not depend on the 
exponent of plasma anisotropy and its equation of state. 



3 Numerical analysis of the energobalance equa- 
tion in Mathematica 

Eg. psp is essentially nonlinear and difficult for analyzing in spite of its appar- 
ent simplicity. Since there is no possibility to find the exact solution of the 
energobalance equation which has important astrophysical applications there is 
a need for its numerical analysis. First attempts of numerical calculation has 
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met significant troubles. Therefore, for the numeric integration control the ana- 
lytic researches were made in [6] . They revealed that the solution has a plateau 
form with minimum value at the point of 7r/2 and after this point the solution 
becomes instable. Also some numerical solutions of the energobalance equation 
in TurboPascal were obtained in the paper. 

Comprehensive analysis of the homogeneous magnetoactive plasma GMSW 
response on a gravitational wave within a wide range of plasma and GW pa- 
rameters was not done in that time due to the software abilities. Nowadays the 
abilities of nonlinear differential equations numeric solution in computer algebra 
system (CAS) Mathematica allow to do such researches. However, the direct 
use of build-in numeric methods towards the energobalance equation ([571) is still 
impossible in case that the governing function derivative changes the sign at the 
point s — tt/2. So one can not change the step size according to the equation 
parameters. 

Empirically was established that the second-order implicit Adams method 
solves the equation much faster and much correct in comparison with other 
explicit and implicit methods. In this research the procedure for numerical 
solving of the differential equation in CAS Mathematica was developed. It 
adapts the integration step according to the parameters T. For all that 
the differential equation is being solving with negative value of derivative up 
to the point 7r/2 using the second-order implicit Adams method. The value 
of function is being taking as the initial value for the positive equation branch 
(|38)) after derivative changes the sign. The integration step is being changed 
and integration method is changing to Euler method which works better in the 
instability region. 

The procedure allows to make analysis of (pS)) numerical solutions depending 
on the first and the second GMSW parameters. It also allows to construct the 
model of magnetoactive plasma response on a GW and to calculate plasma's 
physical characteristics. Numerical researches with our procedure completely 
approve the analytic predictions for the governing function form. At first, the 
solution decreases rapidly then comes to plateau and slowly approaches the 
point of minimum 7r/2 with the function value close to the p9l) . After the point 
of minimum an instability evolves rapidly. For all that the governing function is 
smooth in the whole interval. In Fig[T]the results of numerical solution for the 
energobalance equation in the case of tiny parameter and huge parameter T 
are presented. In this case a process of numeric solution is the most difficult 
and on the other hand the predicted features of the solution are visible. The 
seeming fractures of function by small s and by s = 7r/2 are fake. In fact they 
disappear by scaling up. 

Numeric analysis of the energobalance equation allows to determine the fact 
that the governing function is sufficiently close to the function Ag in the plateau 
area (i.e. the small value of derivative) . Ag nullifies a radical value in right hand 
side of dMl): 




(43) 
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Figure 1. The governing function A(s) by = 10 ^, T = 100, j± = 1/6 

At the point of minimum s — tt/2 this value coincides with the governing 
function minimum value p9p . By increasing the parameter T coincidence of the 
governing function A(s) and the function Ao(s) becomes by the smaller values 
of time s. By small values of the time variable s the governing function is well 
approximated by parabolic law ([55]) . In Fig 15] the plots of the A(s) and the 
Ao(s) functions are shown: 




Figure 2. The A(s) function (sohd line), the Ao(s) function (dashed line), 
the asymptotic if36)) by small values of s: 1 — (dotted line). Everywhere 

= 0.001, T = 10, 7j_ = 1/6. 

This result allows to approximate the plasma response H^/Hq in the plateau 
area of the governing function, i.e. in the maximum response area by the ex- 
pression: 

- — (44) 
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4 Bremsstrahlung response of homogeneous mag- 
net oactive plasma on a gravitational wave 

In a weak GW: 

« 1; L{s) « 1 (45) 

the exact solution of RMHD equations p^ -([T7 | is simplified and we obtain the 
following expressions for physical observed values (see also [S]). 

1. Magnetic field energy density in the comoving frame of reference: 

tt2 tt2 

2. Plasma energy density in the comoving frame of reference: 

e = eoA 1-*=^ . (47) 

3. Physical velocity of plasma: 

l + 2vl l + Ai/27^- ^^^^ 

4. Charged particle density: 

n^noA-i/^'^-L. (49) 

5. Total observed bremsstrahlung intensity detected by resting observer: 

W = WoA-^r^i (^A^/^T-L + A-i/'^'^^) , (50) 
where - total bremsstrahlung intensity in the absence of a PGW [5]: 



Sm^c'^ ° \mc^) ' ^^^^ 

where £ - the kinetic energy of a charged particle. 6. Radiation spectral 
intensity in a high frequencies range where frequencies are comparable with the 
unperturbed cyclotron frequency wj?: 



Q 3ei/o f £o ^ ^ 



2mc V mc 



(52) 



and higher. One may calculate the intensity using standard clectrodynamical 
formulas l8l and find: 



J = JoA-3f(4a-5/2), (53) 
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where 



27r mc^ ' 



Jo 



F{x) = X K5/3{y)dy, 



(54) 



(55) 



K,y{z) - modified Bessel functions of the third kind or Macdonald functions (see 
Ref.[Tn]): 

oo 

ir,(z) = , ^ 

^ ' 2'^r(j/ + i/2) 



e--<^°>^i^t sinh^" t dt. 



In Fig. [3][7] the results of a numeric solution for response of magnetoactive 
plasma on a gravitational wave depending on GMSW parameters are presented. 
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Figure 3. InQuence of the second GMSW parameter T on the relative magnetic 
held energy density evolution H'^/H§ by = 0.0045, 7_l = 1/6.- T = 3 (solid 
line), T = 10 (dashed line), T = 100 (dotted line). 
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Figure 4. InHuence of the second GMSW parameter T on the plasma energy 
density evolution e/eo by ^ 0.0045, 7j_ = 1/6: T = 3 (solid line), T = 10 
(dashed line), T = 100 (dotted line). 
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Figure 5. InRuence of the second GMSW parameter T on the plasma drift 
velocity evolution v^/c by = 0.0045, 7_l = 1/6: T = 3 (solid line), T = 10 
(dashed line), T = 100 (dotted line). 



w 

300r 




Figure 6. InRuence of the second GMSW parameter T on the total observed 
bremsstrahlung intensity evolution W/Wq by — 0.0045, = 1/6: T = 3 
(solid line), T = 10 (dashed line), T = 100 (dotted line). 

Numeric results allow to determine the region of the parameters and T 
where GMSW mechanism becomes sufficiently effective. At that the effect was 
treated as essential if the total observed bremsstrahlung intensity exceeds in 
maximum its initial value in about 2 times. As a result of essential dependence 
of the total observed bremsstrahlung intensity on the parameter this region 
is close to ellipse (Figl?]). 
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Figure 7. Region of the GMSW existence against a quarter ellipse with semi- 
axis 4.2 and 1 background. 
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Figure 8. Bremsstrahlung spectral intensity time evolution in relative units 
by^^ = 0.0045, T = 10 and relative time: s=0; 0.15; 0.31; 0.47; 0.63; 0.78; 0.94; 
1.10; 1.25; 1.41; 1.56^ 7r/2 (from bottom to top). 

In Fig[5] along x-coordinate is a common logarithm of bremsstrahlung fre- 
quency in units of the unperturbed cyclotron frequency wj? and along y-coordinate 
is relative radiation intensity J/ Jo- Maximum of spectral intensity (j53p shifts 
according to the law: 

= 0,29^0 A-5/2(s). 

As was mentioned before an instability rapidly evolves when the governing 
function passes through its minimum which corresponds to the observed values 
maximum. Plasma makes irreversible revers in the direction opposite to the 
PGW propagation direction. This situation is clearly illustrated in FigEl Thus, 
the magnetoactive plasma reacts to a PGW by a single impulse Ref.[2], [3]. 
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Figure 9. Influence of tfie anisotropy parameter on the relative magnetic 
field energy density evolution H'^/Hq by = 0.01, T ~ 10: 7_l = 1/4 (solid 
line), = 1/6 (dashed line), j± = 1/8 (dotted line). 

The numeric analysis results of the anisotropy parameter influence on 
the observed magnetic field energy density are presented in Figl^] One can see 
that resultant anisotropy factor influence is insignificant in spite of the essential 
dependence of the exact solution - on this factor. 

The dependence of the total observed bremsstrahlung intensity semiwidth 
on the GMSW parameters was researched. One can see in FigHU] that the 
GMSW impulse duration i.e. the impulse semiwidth to a high accuracy is equal 
to 7r/4 sa 0.79 or in common units: 

Sr = ^, (56) 

where T is a PGW period. 
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Figure 10. Dependence of the total observed bremsstrahlung intensity semi- 
width W/Wa on T by 7J. = 1/6, = 0.1, = 0.01, ^ o.OOl, = 0.0001 
(from top to bottom). 
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5 Conclusion 



Summarizing the paper results we would like to underline that gravitational 
wave weakness is considered in terms of the conditions (j45[) realization. In this 
case the value of may not be small. Therefore linearity of the theory by 
GW smallness in comparison with 1 (linearity of the Einstein equations left 
hand side by /3) in general does not mean linearity of hydrodynamic theory by 
GW amplitude smallness. Let us notice that such situation is rather unexpected 
though it can be foreseen using a MHD equations exact dimensional analysis. 

The research proved preliminary results of the earlier papers and helped to 
work out in details the GMSW behavior and to describe its evolution process 
in all regions of the parameters. 

Numeric simulation of a magnetoactive plasma response on a gravitational 
wave allows to find next rules of the gravimagnetic shock wave excitation pro- 
cess: 

1. Under realization of the GMSW origin conditions 

C<1; T>1 (57) 

magnetoactive plasma reacts to a PGW by a single impulse where plasma moves 
in the gravitational wave propagation direction. The impulse semiwidth order 
is 1/8 of GW period; 

2. Impulse stops with plasma revers; by this appears the typical impulse form 
(see FigfTTj). It weakly depends on the second GMSW parameter T and is 
determined by the Ao(s) function. 
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Figure 11. Influence of tie first GMSW parameter on the magnetic field energy 
density evolution H'^/H§ by T = 100, 7_l 1/6; = 0.5 (solid line), = 0.3 
(dashed line), = 0.1 (dotted line). 

Under conditions (|57l) the second GMSW parameter T affects only on front 
(small values of time s) and back (time values s are close to 7r/2) edges of the 
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impulse. 

3. Bremsstrahlung spectrum becomes harder during a shock wave passing. 

4. In the maximum of magnetoactive plasma response almost all of a gravita- 
tional wave energy transfers to plasma, magnetic field and to bremsstrahlung 
(see FiginD. 
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Figure 12. Influence of the second GMSW parameter T on tlie GW energy 
absorbtion Aeg/sg = (/3g - /3^)//3o by = 0.0045, jj_ = 1/6; T = 3 (solid 
line), T = 10 (daslied line), T = 100 (dotted line). 
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